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Although most theoretical calculations of quantum wells with non-square profiles assume that 
material composition is varied continuously, it is more common in experiment to grow digital al- 
loys. We compare the Rashba spin-orbit interaction of triangular wells using continuous, discrete, 
and digital alloying profiles in (OOl)-grown triangular InSb/Alf( z )Inii_y>( 2 )Sb, finding a very large 
difference between digital alloying and the others, including a sign change in the Rashba spin-orbit 
coupling. We find that the interface contribution to the Rashba spin-orbit coupling is much larger 
in the continuously- and discretely-alloyed triangular quantum wells than in the digitally-alloyed 
triangular wells, in which it is almost completely absent. The electric field contribution, however, is 
quite similar in all three systems. Due to a much stronger doping dependence in all three systems, 
the electric field contribution dominates at higher dopings, although the very large offset due to the 
near absence of interface contribution in digit ally- alloyed wells persists. 



The spin-orbit interaction couples the spin of an elec- 
tron with the real-space wavefunction, adding a k x in- 
dependent term to the Hamiltonian. This serves as the 
basis for a good deal of very interesting physics, such 
as topological insulators^ and weak antilocalizatiorPEl. 
In semiconductor nanostructures, band structure brings 
about two specific spin-orbit couplings, the Rashba and 
Dresselhaus interactions. Unlike the Dresselhaus spin- 
orbit interaction, which arises from inversion asymmetry 
in the bulk materials, the Rashba interaction arises from 
inversion asymmetry of the structure. This allows the 
Rashba spin-orbit interaction to be controlled via ap- 
plied electric fields, and thus presents a means for co- 
herent control of the electron spin by electrical gates. 
Coherent control of the electron spin through the Ras- 
bha spin-orbit coupling is an important consideration in 
quantum computing with spins in quantum dots^, and 
may also provide the basis for novel devices, such as the 
Datta-Das spin transistor^. 

To grow (001) quantum wells of the form A^( z )Bx_/( z ), 
where f(z) need not be constant, three primary methods 
can be used. The first, continuous alloying, recreates the 
function f(z) by continuously controlling shutter aper- 
ture and cell temperature in the molecular beam epitaxy 
(MBE) system. It is very difficult to tune the parameters 
during growth to reproduce f(z) accurately, however. An 
example of a continuously-alloyed triangular well conduc- 
tion band profile is shown with black solid lines in figure 

m 

It is easier to approximate f(z) by a sequence of steps 
which average f(z) across each step. Growth is then 
stopped between each step, the cell temperatures are al- 
lowed to stabilize, and then growth continues for the step. 
This method, discrete alloying, requires much care, and 
so digital alloying is often used. 

A digitally-alloyed nanostructure takes discrete alloy- 
ing to the next level and builds a nanostructure by creat- 
ing a sequence of square quantum wells The steps are 



then broken into well and barrier regions, with the per- 
centage of well or barrier corresponding to the value of 
f(z) averaged across the step. An example of a digitally- 
alloyed triangular well conduction band profile is shown 
using green dashed lines in figure [T] 

Figure [T] compares continuously and digitally alloyed 
wells and their ground state wavefunctions (from a self- 
consistent calculation described below) in InSb/InAlSb 
and InAs/InGaAs triangular quantum wells. The wave- 
functions are very similar, save for some oscillations in 
the digital alloy. Because of the close similarity in wave- 
functions and greater ease of well growth, digital alloying 
is often used to grow non-constant f(z) wells instead of 
the discrete or continuous alloying methods. 
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FIG. 1: Continuously- and digitally-alloyed well conduction 
band profiles (solid and dashed lines, respectively) and the 
resulting ground states (bottom and top solid lines). Left 
plot is InSb/Alo.2lno.8Sb; right is InAs/Gao^Ino.sAs. The 
first-derivative discontinuities in the InSb wavefunction are 
due to an inhomogeneous effective mass, but do not have a 
large effect on the result; in the InAs well, the effective masses 
are much closer, (color online) 



In this paper, we compare the Rashba spin-orbit in- 
teraction of triangular wells using continuous, discrete, 
and digital alloying profiles. We find that, although the 
electric field component of the Rashba spin-orbit cou- 
pling and ground state wavefunction are similar between 
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continuously-, discretely-, and digitally-alloyed quantum 
wells, the square well basis of the digitally-alloyed well 
causes the interface component of the Rashba coupling 
to be considerably smaller than is the case in the contin- 
uous and discrete wells. 

The Rashba spin-orbit interaction has been showrPESl 
to consist of two parts: an electric field contribution and 
a valence band offset contribution. The electric field com- 
ponent can be understood as an average of the electric 
field by the effective mass wavefunction ip( z ), weighted 
by a bandstructure-determincd factor. This contribution 
takes the form: 



+00 



i^)i 2d #E — 



6+1 



dz b ^ s (E-E Tb {z)-4>(z)Y 



where cJq 



;1 E, 

6m 



(1) 

E Tb (z) is the band edge of the r & 
band through the nanostructure, <f> is the electrostatic 
potential energy from applied and internal self-consistent 
electric fields, ?p(z) is the effective mass wavefunction, E p 
is the k ■ p interband coupling parameter, and mo is the 
bare mass of the electron. 

The other contribution to the Rashba spin-orbit inter- 
action comes from the valence band offsets at the quan- 
tum well interfaces. At a perfectly sharp interface be- 
tween two materials, these offsets form a Dirac delta func- 
tion at the interface. When integrated, this becomes the 
change in the valence bands across the interface and an 
average of the valence band weighting term across the 
interface. This contribution takes the form: 
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where Ni is the number of interfaces, and Zj is the loca- 
tion of the interface i. The interface contribution can be 
understood as the difference between the valence band 
edges at each interface, weighted by the average inverse- 
square of the electrostatic-potential-offset band edge and 
the probability density at the interface. In the case of a 
continuously-alloyed well, the single interface becomes a 
continuous sequence of infinitesimal changes. When used 
on a real-space grid, as in this work, the material is con- 
stant within each gridsite. This yields a discrete alloy 
with steps of one gridsite. 

Examining a square well with a constant electric field 
(fig. [2]) can provide intuition. In this square well, the 
valence band difference at the left and right interfaces 
are of equal magnitude but opposite sign. The inverse- 
square band offset weighting terms are also the same on 
the left and the right, if we neglect the electric field term. 
If a non-uniform electrostatic potential were introduced, 
for instance from the internal electric field between elec- 
trons in the well and the dopant layer, a small difference 
would arise. A much larger contribution comes from the 
difference in probability density at the left and right in- 



Energy 
U~ 2 



FIG. 2: Cartoon of a wavefunction in the conduction band 
of a square well. A constant electrostatic field 4>(z) = kz has 
been applied to the structure. Valence bands are not shown 
due to space considerations, but are also square. 



terfaces. This will in general be the largest contribution 
to the interface term. 

It should be noted that these equations do not account 
for Dirac renormalizatiorfiil and thus neglect e.g. the 
Darwin contribution to the overall Hamiltonian. This is 
not a major impediment to the current discussion, as the 
Dirac renormalization terms have only a modest effect on 
the above equations. 

Much more interesting is the energy E above. The 
effective mass wavefunction perturbation methods used 
to determine equations [T] and [2] are derived using a pro- 
jection operator formalism which has been described by 
LowdinP^l Ultimately, the E term in equations [l] and [2] 
comes from the solution to the full N-band Schrodinger 
equation for all of the carriers in the well. Both aspects of 
this are important to fully understand the Rasbha spin- 
orbit interaction in doped quantum wells: the fact that 
these come from a solution to the full Schrodinger equa- 
tion as well as the fact that it must include every carrier 
in the well. 

One approach is to use the energy of the particular 
subband being examined. In this approach, the single- 
particle, effective mass ground state eigenenergy (Eq) is 
used for E. The drawback of this approach is that it 
neglects the carriers at higher k, at energies up to the 
Fermi energy. This approach is taken by some groups^, 
although they approximate Eq by the energy of the con- 
duction band edge in the well. Although Eq is readily 
understood in triangular or more general quantum wells, 
the value to take for the conduction band edge in trian- 
gular or more complicated wells is unclear. 

Alternately, a clear second approach is to take the 
Fermi energy {Ep) for E. This then can be qualita- 
tively understood as examining only the highest occu- 
pied state. This approach has also been taken by a num- 
ber of groups^Hmj, anc i like the first approach, som^ 
approximate Ep with the conduction band edge in the 
well. However, other groups use Ep. In addition, using 
the Fermi energy instead of the conduction band edge in 
the welP should qualitatively reproduce the trend of de- 
creasing Rashba spin-orbit coupling with increasing car- 
rier density (i.e. increasing Ep) due to the inverse-i? 2 
dependency of a, which is harder to explain with other 
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means. 

In this paper, we use the latter definition of E, namely 
that we simply replace E with Ep. However, simply in- 
troducing Ep is somewhat ad hoc and a more rigorous 
formalism is certainly called for to clarify the ambigu- 
ity in methods. Koga et alP compared the approxima- 
tion for E used by Engels et alP^l (the conduction band 
edge in the well) and found approximately a 20-30% in- 
crease in a when using Ep instead of the band edge in 
InAlAs/InGaAs/InAlAs wells (70-90 meV). The ground 
state energy, necessarily being lower than Ep and higher 
than Eq should lie somewhere below that. The difference 
in InSb/InAlSb is likely even larger, due to the smaller 
bandgap in InSb. Using the Eq or conduction band edge 
approximation of E should change only the scale of our 
results, not the results themselves. 

To compare digitally-, discretely-, and continuously- 
alloyed triangular wells, we have calculated the 
Rashba spin-orbit coupling in doped continuously- 
alloyed, discretely-alloyed, and digitally-alloyed quantum 
wells without an applied electric field. The wavefunctions 
within the conduction band were calculated using a real- 
space self-consistent effective mass method, using an in- 
homogeneous, material-dependent effective mass. As the 
abrupt change in effective mass between the well and bar- 
rier in InSb/AlInSb quantum wells causes discontinuities 
in the first derivative of the wavefunction, we additionally 
calculated the same well types but with an InAs/GalnAs 
quantum well system where the effective masses in the 
well and barrier are much closer. 

To prevent infinite loops in the self-consistent calcula- 
tion, both the second-to-last and previous wavefunctions 
were compared to the current wavefunction to determine 
convergence, with no potential mixing being performed. 
An alternate method would have been to mix a fraction 
of the previous and current potential values.^ 

After the effective mass wavefunctions were calculated, 
equations[l]and[2]were used to calculate the Rashba spin- 
orbit coupling for the conduction electrons in the ground 
state subband. 

The Rashba couplings at a variety of dopings were 
thus calculated and are plotted in figure [5) The elec- 
tric field and interface contributions are plotted in fig- 
ures [3] and [4] respectively. Doping was varied be- 
tween 2 x 10 cm~ 3 and 1.2 x 10 18 cm~ 3 across 2.5nm 
in the InSb/Iny^Al^^^Sb quantum wells, and be- 
tween 4 x 10 16 cm -3 and 2 x 10 17 cm -3 across 2.5nm in 
the InAs/Iny( z )Gai_/( z ) As quantum wells. The doping 
regions in both cases began 15nm to the right of the 
sharp interface of the triangular wells. The wells were 
all 28.8nm wide, using a grid spacing of 0.4nm for the 
continuously and discretely alloyed wells and O.lnm for 
the digitally alloyed well. 

Due to the similarity in wavefunction shape, the elec- 
tric field term can be expected to be very similar for 
digitally- and continuously-alloyed wells. This expecta- 
tion is borne out in our calculations, shown in subfigures 
3(a)| and [3(b)] There is a strong increase in the mag- 
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FIG. 3: Electric field contribution to the Rashba 
spin-orbit coupling in InSb/Alo.2lnrj.8Sb (left plot) and 
InAs/Ino.2Gao.8As (right plot) continuously-alloyed (solid 
line), discretely-alloyed (dotted line), and digitally-alloyed 
(dashed line) triangular quantum wells. Results for discrete 
and continuous alloys are nearly indistinguishable. 



nitudc of the Rashba spin-orbit interaction as doping is 
increased for all three alloying types. This is due to the 
increasing electric field between the carriers in the well 
and the dopant layer as doping is increased, and is ex- 
pected. 
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FIG. 4: Interface contribution to the Rashba spin-orbit cou- 
pling in InSb/Alo.2lno.8Sb (left plot) and InAs/Ino.2Gao.8As 
(right plot) continuously-alloyed (solid line), discretely- 
alloyed (dotted line), and digitally-alloyed (dashed line) tri- 
angular quantum wells. 



In a continuously alloyed well, given that one side is 
a sharp interface, the other side of the well provides a 
continuous (or monotonically changing) set of interfaces. 
In a digitally alloyed well, by contrast, the sequence of 
square wells present both a strong positive and negative 
interface at the subwell boundaries. This can be seen 
in figure [T] where the continuously-alloyed and digitally- 
alloyed wells and wavefunctions are superimposed (with- 
out the internal potential) for InSb and InAs triangular 
quantum wells. Thus instead of one large change in the 
valence band fighting a large set of small changes, a dig- 
itally alloyed well has two large changes in opposition. 
In addition, figure [T] shows an increasing probability den- 
sity in the first subwell (going right to left) which is then 
countered by a series of subwells with decreasing prob- 
ability density. Because the valence band changes are 
identical in each of the subwells, this increase in proba- 
bility density in the first subwell followed by a decrease 
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across the subsequent subwells seen in figure [T] serves to 
counter itself in addition to the small changes within each 
individual subwell. 
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FIG. 5: Total Rashba spin-orbit coupling in 

InSb/Alo.2lno.8Sb (left plot) and InAs/Ino.2Gao.8As (right 
plot) continuously-alloyed (solid line), discretely-alloyed 
(dotted line), and digitally-alloyed (dashed line) triangular 
quantum wells. Continuous and discrete alloys are almost 
indistinguishable. 



There is a moderate increase in the magnitude of the 
interface contribution to the Rashba spin-orbit interac- 
tion as doping is increased. This has its origin in two 
parts: the first is the increase in electrostatic potential 
on one side of the well due to the internal electric field 
between the carriers in the well and the dopant layer, 
leading to a larger asymmetry across the wells in the 
valence band. The second part is the wavefunction tun- 
neling into the barrier increasing as the internal electric 
field increases, leading to another increase in the asym- 
metry of the well and therefore an increasing interfacial 
Rashba contribution. 

The sum of the two contributions is plotted for the 
quantum wells in figure [5j Because of the almost com- 
plete absence of the interface contribution in the digital 
alloys but the similarity of the electric field contribution 



to that found in discrete or continuous quantum wells, 
the total Rashba spin-orbit interaction in the digitally 
alloyed wells is considerably different from the continu- 
ous and discrete counterparts. This is perhaps surprising, 
due to the similarity of digital and continuous alloy wave- 
functions. This further reinforces the importance of the 
interface term in Rashba calculations. 

Aside from the offset due to the differences in inter- 
face terms, the similarity in electric field terms in digital 
compared to discrete and continuous alloying causes the 
high-doping regimes, where the electric field term domi- 
nates, to have similar trends. This has its origin in the 
greater sensitivity of the electric field contribution to the 
internal electric field compared to the interface contribu- 
tion as doping is increased. 

To summarize, digitally alloyed quantum wells have 
substantially different Rashba spin-orbit couplings com- 
pared to continuously or discretely alloyed wells, even 
changing sign in some doping regimes. This has its origin 
in the near absence of the interface contribution in digital 
alloys. As the electric field contribution is more sensitive 
to doping than the interface contribution, the trend as 
a function of doping density in the high-doping regime 
in digital alloys are similar to the trend in continuously 
or discretely alloyed wells. This is likely to be an im- 
portant consideration in further experimental growth of 
novel quantum wells for topological insulators and other 
applications. 
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